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Abstract 

Several experiments measure the fluorescence light produced by extensive air show- 
ers in the atmosphere. This light is converted into a longitudinal shower profile 
from which information on the primary energy and composition is derived. The 
fluorescence yield, as the conversion factor between light profile measured by EAS 
experiments and physical interpretation of showers, has been measured in several 
laboratory experiments. The results, however, differ considerably. In this article, a 
model calculation of the fluorescence emission from relevant band systems of nitro- 
gen in dependence on wavelength and atmospheric conditions is presented. Different 
calculations are compared to each other in combination with varying input parame- 
ters. The predictions are compared with measurements and the altitude-dependence 
of the fluorescence yield is discussed in detail. 
PACS: 96.40.Pq 

Key words: fluorescence light emission; extensive air shower; atmosphere 



1 Introduction 



This is the second article of a series of investigations of the importance of atmo- 
spheric properties for the reconstruction of extensive air showers (EAS). The 
first article [1] describes the effect of changing atmospheric density profiles on 
the longitudinal EAS development. In particular, if the fluorescence technique 
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is applied, the conversion of atmospheric depth, as used in shower simulations, 
to geometrical altitude, as reconstructed from fluorescence measurements, and 
vice versa is very important for primary cosmic ray mass reconstruction. 

This article addresses the fluorescence light emission of EAS which is used 
for the determination of the total energy of EAS. In several air shower ex- 
periments, for example, HiRes [2], the Pierre Auger Observatory [3,4,5], and 
Telescope Array [6,7], the fluorescence technique is employed for detecting 
EAS. Measuring the fluorescence light that nitrogen molecules emit after be- 
ing excited by charged particles traversing the atmosphere is the most direct 
method of detecting the longitudinal shower profile. For the event reconstruc- 
tion procedures of these air shower experiments, the knowledge of the fluo- 
rescence yield FYx and its dependence on atmospheric conditions are crucial 
parameters. 

The Pierre Auger Observatory is up to now the only existing EAS experi- 
ment which applies hybrid detection techniques. The secondary particles of an 
EAS are measured at ground and simultaneously the fluorescence light of the 
longitudinal shower development is detected with telescopes. For extracting 
a cosmic ray spectrum from the data, the events detected with ground de- 
tectors are analyzed while the energy calibration is deduced from fluorescence 
detector events and the correlation of these two types of events is derived from 
hybrid events [8]. This cosmic ray spectrum and its comparison with spectra 
published by other experiments (AGASA [9] and HiRes [10], for a recent re- 
view see [11]) reveal that the fluorescence yield might be a crucial parameter 
for the energy reconstruction of air showers. For the conversion of detected 
fluorescence light to energy deposited in the atmosphere by EAS and flnally 
to the total energy of the primary particle, not only the total fluorescence 
yield in the detected wavelength (A) region is important but also the spectral 
distribution. For example, the emitted light suffers Rayleigh scattering while 
traversing the atmosphere towards the telescopes. Since the scattering cross 
section has a A^^ dependence, the long-wavelength part of the fluorescence 
spectrum has a higher transmission than the short-wavelength region. 

The outline of this paper is as follows. In Section 2, the fluorescence emission 
in air is reviewed and an analytical model (Sec. 2.1) for calculating the flu- 
orescence emission in dependence on wavelength and atmospheric conditions 
is described. Particularly, the band systems of nitrogen contributing mainly 
to the fluorescence light emission are discussed. A compilation of several pa- 
rameters from different authors used in calculations of this paper is given in 
Sec. 2.2. In Sec. 3 a detailed comparison of measurements with the calculations 
is done. The aim is to combine the laboratory measurements with our current 
understanding of the fluorescence emission processes in the atmosphere, to 
show possible sources of uncertainty, and to provide an easy way of account- 
ing for varying atmospheric conditions. The dependence on these atmospheric 
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conditions is explicitly presented in Sec. 4. 



2 Model for Fluorescence Light Emission in Air 

The most numerous charged particles in an EAS are electrons and positrons. 
Their energy deposit in air by ionization and excitation of air molecules gives 
rise to fluorescence light emission. In the wavelength region between 300 and 
400 nm, the band systems with the strongest emissions are found. Therefore, 
all EAS experiments using the fluorescence technique apply UV fllters with 
largest transmittance between roughly 310 and 400 nm. The residual wave- 
lengths are cut in favor of reducing the night sky background. 

The major components of the atmosphere are 78.08% N2, 20.95% O2, and 
0.93% Ar per volume. All three constituent parts influence the emission of 
fluorescence light, however, with strongly differing importance. 

The main fluorescence light is emitted by two electronic states of N2, these 
are the second-positive (2P) band system, C'^Uu - -B^H^, and the first-negative 
(IN) system of Ng", B'^^t ~ -^'^^t ' ^^'^^ band within a system belongs to a 
transition from a vibrational level of the upper state i>' to a vibrational level 
of the lower state u" [18]. The band structure is caused by the rotational 
levels and consists of a large number of spectral lines which are very close to 
each other [19]. These lines will not be distinguished in this work. The bands 
usually have a band head at one end where the intensity falls off suddenly. In 
the following, the wavelength given for each band belongs to the position of 
the corresponding band head. 

For both, absorption and emission, the Franck- Condon principle applies in 
good approximation. Three processes of excitation of N2 can be distinguished: 

• Direct excitation: The energy deposited in air excites nitrogen molecules 
proportional to an energy-dependent cross section a^'{E). This process 
mainly acts on the N2' IN system 



• Excitation via secondary electrons: High energy particles in EAS ionize N2 
producing several lower energy secondary electrons. These e~ are able to 
excite also the N2 2P system with a resultant spin change 



N2 + e ^ N+* + e + e. 



(1) 



N2 + e(T)^N;(C3nj + e(i). 



(2) 



In addition, the 2P system can also be excited by recombination 



N+ + e ^ N;(C' 




(3) 



3 



• Via Auger electrons: Some ionization processes will release K-electrons lead- 
ing to the emission of Auger electrons. These are on their part again able 
to excite the N2 molecules. However, it must be emphasized that the cross 
section for this K-shell ionization is much lower than those for the processes 
mentioned above. The cross sections for N2 excitation and for the ionization 
to N2' are of the order of 10~^^ m^ [20,21] and for the K-shell ionization of 
the order of lO'^^ [22]. 

In air, the optical emission of the prompt radiative return from the upper 
states of the 2P and IN system of nitrogen will be affected by some compet- 
ing processes. The most important process is coUisional quenching. Excited 
nitrogen molecules might collide with other molecules in air before the de- 
excitation via fluorescence light emission happens. 

Argon can be excited by the reaction e-|- Ar — > Ar*, where the excitation cross 
section is largest for Ar(^P2) [12]. This process is followed by Ar* + N2 ^ 
Ar + N2(C^n„). The energy is mainly transferred from argon to nitrogen via 
secondary electrons rather than direct collisions [13]. The excited state C^Hu 
is the upper level of the second-positive (2P) system of N2 which radiates 
photons mainly in the wavelength region between 300 and 400 nm. The lower 
state is B^Ug. This increase of the emission competes, however, with a higher 
quenching rate, which leads to non-radiative de-excitation, due to additional 
collision partners in the form of argon atoms in air. The net effect of argon 
is expected to be less than 1% contribution to the fluorescence hght [14]. 
However, argon emits also directly fluorescence light at around 310 nm [15]. 
This transition, A'^T,+ - X^H, has been investigated in argon water- vapor 
mixtures and the highest intensity has been found for very low argon pressure 
and 0.06 hPa water vapor. For EAS experiments, this contributions will be of 
minor importance, too. 

The UV- fluorescence light emission from O2 is negligible [16]. The contribution 
in the relevant wavelength region stems from O^A^H^ — X^Ug transmissions. 
However, already the Einstein coefficients ^ are reduced in average by a factor 
of about 30 compared to the emissions of the 2P system of N2 [17]. The 
emissions of atomic oxygen start at 395 nm and go up to 845 nm [16]. These 
bands are of no importance for EAS experiments. 

Generally, it is assumed that the fluorescence light is proportional to the energy 
deposit of an EAS. The contribution of electrons and positrons to the energy 
deposit according to the initial kinetic energy distribution in an air shower 
has been studied elsewhere [23]. Only 10% of the energy deposit stems from 
particles with energies less than 0.1 MeV. Particles with energies between 
0.1 and 10 MeV contribute 35%, between 10 and 100 MeV also 35%, and 



^ radiative transition probabilities 
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between 100 and 1000 MeV 17%. The remaining 3% originate from particles of 
energy above 1000 MeV. Several experiments have recently started to measure 
the proportionality of fluorescence yield to energy deposit. So far none of 
the results on the energy dependence is absolutely cahbrated, thus a direct 
comparison is difficult. First relative results seem to confirm the expected 
correlation in different energy regions, see e.g. at about 1 MeV Nagano et 
al. [24], between 50 and 420 MeV Bohacova et al. [25], and at 28.5 GeV Belz 
et al. [26]. 

Depending on their initial energy, EAS particles produce secondary electrons 
of various lower energies. These can excite the N2 but they may also suffer 
an attachment process: if, on their way from the production site to the N2 
molecules, the secondary electrons encounter a strong electronegative pollu- 
tant (e.g. O2, H2O, CO2, H2, Xe, CH4 which are trace gases in the atmosphere), 
they are attached to this pollutants and cannot excite the N2 molecules any- 
more [13]. This process is beyond the scope of this paper. 

2.1 Mathematical Description 

The existing results of fluorescence yield measurements show quite large dif- 
ferences. Furthermore, the data have to be applied to air shower reconstruc- 
tion procedures. While secondary particles of EAS traverse from high to low 
altitudes, they encounter continuously changing atmospheric conditions. Ad- 
ditionally, the atmospheric conditions vary from day to day with the largest 
differences between the seasons summer and winter at the sites of all ex- 
isting air shower experiments. The aim of the calculations shown here is to 
cross-check the laboratory measurements with the current understanding of 
the processes in the atmosphere, show possible sources of uncertainties, and 
provide an easy way to account for varying atmospheric conditions in EAS 
measurements. 

The quantum efficiency of fluorescence can be defined as 

rate of de-excitation via radiation Tc , . 



where the rate of de-excitation is proportional to the reciprocal of the life 
time. The mean life time of the radiative transition to any lower state is tq 
and to collisional quenching Tc. The collisional quenching can be described by 
kinetic gas theory. The molecules, in the case of air, move with velocities fol- 
lowing the Maxwell-Boltzmann distribution which is strongly correlated with 
gas temperature. As a good approximation, the collision rate depends on the 
mean velocity of molecules v — The resulting mean hfe time due to col- 



total rate of de-excitation 



To + Tc ' 



(4) 
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lisional quenching is the ratio of the mean free path, in this case for molecules 
of one type moving with roughly the same velocity, and the mean velocity: 



Te = (V2 • P„ • (7NN • V)-^ = ^J{7^M/kT) ■ (4p„ • (7nn)~', (5) 

where pn is the particle number density, ctnn the coUisional cross section be- 
tween nitrogen molecules, T the temperature, k the Boltzmann constant, and 
M the molecular mass. Now the pressure dependent fluorescence efficiency can 
be written as 



with £^ being the fluorescence efficiency at wavelength A without coUisional 
quenching, n denoting the number of photons, the energy of a single photon 
with the corresponding wavelength, E^ep the deposited energy in the observed 
medium, and p/p^i = r^y /T^y. The pressure p is that of the observed medium 
(e.g. air), p^, is a reference pressure at which tq is equal to Tc- Toy and Tcy 
are the mean life times for excitation level u'. Applying actual atmospheric 
conditions, with air taken to be a two-component gas, the relation between p 
and p'j^, can be written as 



^ = '^oy ■ ( J r H r ) (7) 

Pu' VTNN,i/'(,C"NN,i/'j TnO,i/'(,CNO,i/'J- 



To,.'Pair • Na kTNA 



R-T 



-^■U-C.{^2)-a^^yJ-— (8) 

TT V V ^^^m,N 



with Avogadro's number Na, the masses per mole for nitrogen Mm,-H and 
oxygen M„i^Q, the universal gas constant R, the cross sections for coUisional 
de-excitation for nitrogen-nitrogen a^-^y and nitrogen-oxygen (T-^^oy, and the 
fractional part per volume Cy of the two gas components. 



2.2 Input Parameters 



To estimate the fluorescence efficiency with these equations, several parame- 
ters have to be obtained from measurements and/or calculations. Most impor- 
tant is the fluorescence efficiency without coUisional quenching e\. An early 
compilation of elder measurements performed by Bunner [14] provides these 
values for 18 band systems of the 2P and for 1 band system of the IN state of 
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nitrogen in the wavelength region between 300 and 400 nm. In a more recent 
pubhcation by Gilmore et al. [17], the Einstein coefficients Aj,/,^// of the 2P 
nitrogen state for the transitions from i/' = ... 4 to i/" = ... 21 and the ra- 
diative life times Toy ior u' — . . A are given. For the IN nitrogen state, the 
Einstein coefficients A,^i^'/ for the transitions from i/' = ... 10 to i/" = ... 21 
and the radiative life times Toy for u' — ... 10 are listed. The intensity of 
a transition could be calculated by N^/^i^n = Toy ■ A^^ij^n ■ N*,, where N*, is 
the number of excited states. Since this number is unknown, a relative fluo- 
rescence efficiency is calculated by multiplying the Einstein coefficients with 
the radiative life times and a relative apparent excitation cross section Qapp 
[27] . The apparent excitation cross sections are the sum of the optical emission 
cross sections over v" for a given v' . The excitation cross sections are called 
apparent, because they represent the sum of direct excitation and cascades 
from higher levels down to that particular one. The Qapp values are derived 
from the apparent excitation cross sections by normalizing the value for u' — 
to unity of the electronic state. Values are given in the publication by Pons 
et al. [27] for the 2P band system of nitrogen for u' = ... 4 and for the IN 
band system for i/' = ... 3 in Stanton and St. John [28]. 

The relative fluorescence efficiency can then be normalized for each electronic 
state to, e. g., the most prominent band. For the 2P system of nitrogen, it is 
the (0-0) band with a wavelength of 337.1 nm, and for the IN system it is 
again the (0-0) band with the corresponding wavelength of 391.4 nm. In our 
calculations, the efficiency values for these wavelengths as given in Bunner [14] 
are chosen. 

In the following the fluorescence efficiencies are calculated for the bands of the 
2P system. They dominate the fluorescence light emission in the atmosphere 
due to excitations caused by EAS. The fluorescence efficiencies labeled with 
Sx,G.-F. are obtained by applying the Einstein coefficients given by Gilmore 
et al. and the relative apparent excitation cross sections from Eons et al. The 
values used within this paper are listed in Tab. 1. 

Further parameters in the upper equations are the deactivation constants 
which are the radiative life time Toy and the coUisional cross section be- 
tween nitrogen and nitrogen molecules cr^Ny and between nitrogen and oxy- 
gen molecules (TNoy- The values for Toy obtained by Gilmore et al. are listed 
in Tab. 1 and 3. Bunner provides collisional cross sections and radiative life 
times for the most prominent band systems of nitrogen, see Tab. 2. Recent 
measurements by Morozov et al. [29] were performed for the 2P u' — 0,1 band 
systems, see also Tab. 2. In further calculations presented in this article, la- 
beled with Morozov, the values from Bunner are replaced by the newer data 
by Morozov et al. where available. 

No coUisional cross sections are available for the first-negative system of nitro- 
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Table 1 



Constants for the fluorescence eflficiency £° of the 2P system of nitrogen. 



Band 


Wave- 




Gilmore et al. [17] 


Pons 








length 


\,Bunner 




Toy 


et al. [27] 


^\,G.-F. 


{u'- 


-v") 


A (nm) 


(%) 




(s) 


Qapp 


(%) 


2P 


(0-0) 


337.1 


.082 


1.31E7 


3.71E-8 


1 


.082 


2P 


(0-1) 


357.7 


.0615 


8.84E6 


3.71E-8 


1 


.0553 


2P 


(0-2) 


380.5 


.0213 


3.56E6 


3.71E-8 


1 


.0223 


2P 


(1-0) 


315.9 


.050 


1.19E7 


3.75E-8 


0.7 


.0527 


2P 


(1-1) 


333.9 


.0041 


5.87E5 


3.75E-8 


0.7 


.0026 


2P 


(1-2) 


353.7 


.029 


5.54E6 


3.75E-8 


0.7 


.0245 


2P 


(1-3) 


375.5 


.0271 


4.93E6 


3.75E-8 


0.7 


.0218 


2P 


(1-4) 


399.8 


.016 


2.43E6 


3.75E-8 


0.7 


.0108 


2P 


(2-1) 


313.6 


.029 


1.01E7 


3.81E-8 


0.26 


.0169 


2P 


(2-2) 


330.9 


.002 


7.99E5 


3.81E-8 


0.26 


.0013 


2P 


(2-3) 


350.0 


.004 


1.71E6 


3.81E-8 


0.26 


.0029 


2P 


(2-4) 


371.0 


.010 


4.04E6 


3.81E-8 


0.26 


.0068 


2P 


(2-5) 


394.3 


.0064 


3.14E6 


3.81E-8 


0.26 


.0052 


2P 


(3-2) 


311.7 


.005 


5.94E6 


3.90E-8 


0.081 


.0032 


2P 


(3-3) 


328.5 


.0154 


2.85E6 


3.90E-8 


0.081 


.0015 


2P 


(3-4) 


346.9 


.0063 


1.15E5 


3.90E-8 


0.081 


.0001 


2P 


(3-5) 


367.2 


.0046 


2.35E6 


3.90E-8 


0.081 


.0013 


2P 


(3-6) 


389.5 


.003 


3.00E6 


3.90E-8 


0.081 


.0016 


2P 


(4-3) 


310.4 




3.02E6 


4.04E-8 


0.041 


.0008 


2P 


(4-4) 


326.8 




3.71E6 


4.04E-8 


0.041 


.0010 


2P 


(4-5) 


344.6 




1.24E5 


4.04E-8 


0.041 


.0000 


2P 


(4-6) 


364.2 




9.98E5 


4.04E-8 


0.041 


.0003 


2P 


(4-7) 


385.8 




2.33E6 


4.04E-8 


0.041 


.0007 
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Table 2 



Deactivation constants for air in the lower atmosphere. 





Bunner [14] 




Morozov et al. 


29] 




0"NO 


<7NN 


tq 


<7NN 




To 




(m2) 


(m2) 


(s) 


(m2) 


(m2) 


(s) 


IN = 


13E-19 


4.37E-19 


6.58E-8 








2P z/ = 


2.1E-19 


l.OE-20 


4.45E-8 


1.82E-20 


8.53E-19 


4.17E-8 


V = 1 


5.0E-19 ^ 


3.5E-20 


4.93E-8 


3.77E-20 


8.04E-19 


4.17E-8 


V = 2 


7.0E-19 ^ 


8.8E-20 


4.45E-8 








u = 2, 


8.0E-19 


1.2E-19 


6.65E-8 









^ This value is determined by the given results of [14] and not given in the original 
publication. 



gen for higher excitation levels than the (0-0) band. For an estimation of higher 
level contributions, one can assume that the coUisional cross section between 
the ionized nitrogen molecule and further nitrogen and oxygen molecules are 
equal to that of the (0-0) transition. Since the coUisional cross section usually 
increases for higher excitations, this is an assumption which gives an upper 
bound on the fluorescence light from these levels. The relative fluorescence 
yields for the IN bands arc estimated with the same method as already used 
for 2P, however, applying the parameters given in Gilmore et al. and Stanton 
and St. John (see e\ Q.-st. 3). Even under these conditions, the IN (1-0) 

and IN (1-1) transitions add up to less than 1% to the total fluorescence yield 
between 300 and 400 nm at sea level in the US Standard Atmosphere. Higher 
excitations give even smaller contributions. Therefore, we will only consider 
the (0-0) transition of the IN band and neglect the higher excitations. 

The calculation introduced in Sec. 2.1 together with the ^A^Mnner ^^e 
deactivation constants from Morozov is used as reference model in this paper. 
It benefits from the completeness of the Bunner data and from the accuracy 
of the measurements from Morozov et al. 



3 Compcirison with Measurements 

Wavelength-dependent results of fluorescence yield measurements have been 
provided in three publications [14,30,31]. Bunner lists several intermediate 
values: £x^'{p^ T) in %, and the fluorescence efficiency e^^' in units of pho- 
tons/MeVof deposited energy which is ef''-{p,T)-{\/hc) at sea level (s./.). The 
values for £^'' (p, T) and SE^ep gi'^^n explicitly in [14] are not reproduced by the 
calculations shown here, see Table 4. Possible reasons are rounding uncertain- 
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Table 3 



Constants for the fluorescence efficiency of the IN system of nitrogen. 



Band 


Wave- 




Gilmore et al. [17] 


Stanton &: 






length 


X,Bunner 




Toy 


St. John [28] 


^X,G.-St. 


{u' - v") 


A (nm) 


(%) 


(1/s) 


(s) 


Qapp 


(%) 


IN (0-0) 


391.4 


.33 


1.1 1E7 


G.23E-8 


1 


0.33 


IN (1-0) 


358.0 


- 


5.76E6 


6.20E-8 


0.117 


0.0193 


IN (1-1) 


388.2 


- 


4.03E6 


6.20E-8 


0.117 


0.0135 


IN (2-0) 


330.5 




9.02E5 


6.19E-8 


0.009 


0.0002 


IN (2-1) 


356.1 




7.88E6 


6.19E-8 


0.009 


0.0021 


IN (2-2) 


385.5 




9.27E5 


6.19E-8 


0.009 


0.0003 


IN (3-1) 


329.6 




2.08E6 


6.23E-8 


0.004 


0.0002 


IN (3-2) 


354.6 




8.09E6 


6.23E-8 


0.004 


0.0009 



ties by Bunner or the use of deviating numbers for variables concerning air 
conditions. Davidson and O'Neil [31] list results for £^' (]9, T) for wavelengths 
above 320 nm. It should be mentioned that the results in [31] are given for 
p — 800 hPa. The increase of the total fluorescence yield between 300 and 
400 nm from sea level with p = 1013 hPa to approximately 2 km a.s.l. with 
p = 800 hPa amounts to about 2%. Nagano et al. report directly the values 
for FYx at sea level for 0.85 MeV electrons [30], however, only 10 contributing 
emission bands are listed. For comparing the results of all authors, 0.85 MeV 
electrons are chosen as exciting particles, so the ionization energy deposit is 
dE/dX = 0.1677 MeV/kg-m~^ [24]. It is assumed that the fluorescence yield 
is proportional to the energy deposit as discussed in Sec. 2. Air is taken to 
be a composition of 78.8% N2 and 21.1% O2 per volume [24]. The resulting 
fluorescence yield can be written as 



Tpv / >^ dE rphotonsl 

FYx = T)- — ■ — ■ Pair ■ 9) 

he dX L m J 

A comparison of the obtained FYx values at sea level in the US Standard 
Atmosphere (US-StdA) [32,1] is shown in Table 4 and illustrated in Fig. 1. 

The total fluorescence yield reported by Bunner directly in [14] is much lower 
than the other measurements and calculations. The total value from David- 
son and O'Neil is higher by 6.8% for wavelengths above 320 nm as compared 
to our model for the same wavelength region. The calculations shown here 
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Table 4 



Fluorescence yield at sea level in the US Standard Atmosphere. For comparing the 
results of all authors, 0.85 MeV electrons are chosen as exciting particles, so the 
ionization energy deposit is 0.1677 MeV/kg-m~^. See text for details. 





PY^.l. ^ photons^ measurements from 


FY^^- (Photons-) calculations with 


Wave- 


Bunner 


Davidson & 


Nagano 


\, Bunner ' 


\, Bunner ' 




length 


[14] 


O'Neil [31] 


et al. [30] 


Table 2 - 


Table 2 - 


'^Nx, Morozov: 


A (nm) 








Bunner 


Morozov 


TO,Gilmore 


310.4 


_ ^ 


_b 


- ^ 


_ 


_ 


0.001 


311.7 


0.008 


_b 


_ a 


0.009 


0.009 


0.010 


313.6 


0.090 


_b 


_ a 


0.094 


0.094 


0.064 


315.9 


0.224 


_b 


0.549 


0.240 


0.279 


0.326 


326.8 


_ ^ 


_ ^ 


_ a 


_ 




0.002 


328.5 


0.027 


0.035 


0.180 


0.029 


0.029 


0.005 


330.9 


0.007 


_ a 


_ ^ 


0.007 


0.007 


0.005 


333.9 


0.019 


_ a 


_ a 


0.021 


0.024 


0.017 


337.1 


0.887 


1.173 


1.021 


1.169 


1.108 


1.242 


344.6 


_ a 


_ a 


_ a 


_ 


_ 


0.000 


346.9 


0.012 


0.015 


_ a 


0.012 


0.012 


0.000 


350.0 


0.014 


0.013 


_ a 


0.014 


0.014 


0.012 


353.7 


0.146 


0.188 


0.130 


0.156 


0.181 


0.170 


357.7 


0.707 


0.889 


0.799 


0.931 


0.882 


0.889 


364.2 


a 


a 


_ ^ 


_ 


_ 


0.001 


367.2 


0.009 


0.012 


_ a 


0.010 


0.010 


0.005 


371.0 


0.037 


0.047 


_ a 


0.038 


0.038 


0.030 


375.5 


0.150 


0.187 


0.238 


0.155 


0.180 


0.160 


380.5 


0.261 


0.328 


0.287 


0.343 


0.325 


0.381 


385.8 


_ a 


_ a 


_ a 






0.001 


389.4 


0.006 


_ a 


_ ^ 


0.007 


0.007 


0.006 


391.4 


0.281 


0.454 


0.302 


0.315 


0.315 




394.3 


0.025 


0.032 


0.063 


0.026 


0.026 


0.025 


399.8 


0.090 


0.119 


0.129 


0.097 


0.113 


0.085 



This transition has not been measured. 
Only measurements above 320 nm. 
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continuation Table 4 


Wave- 
length 
A (nm) 


pys.l. (I 

Bunner 
[14] 


)hotons \ jj-^gjj^g^-^j 
m ' 

Davidson & 
O'Neil [31] 


cmcnts from 
Nagano 
et al. [30] 


^0 

\, Bunner ' 

Table 

9 

^Bunner 


lotons \ 1 

calcu 

111 ' 
\, Bunner ' 

Table 

•^M orozov 


lations with 

^A,G.-F.' 
^Nx,Morozov i 
TOyGilmore 


sum of ^ 


= (300-4 
3.001 


30) nm 
3.490 ° 


3.698 


3.672 


3.653 


3.438 ^ 


sum of a 


11 Nagano 
2.798 


-wavelengths 
3.405 


3.698 


3.460 


3.438 


3.283 


sum of a 


11 Nagano 
2.574 


-wavelengths a 
3.405 


bove 320 nm 
3.149 


3.221 


3.159 


2.957^ 


sum of a 


11 Nagano 

2.293 


-wavelengths a 

2.951 


bove 320 nm { 

2.8:1.7 


md without 

2.906 


391.4 nm 

2.811 


2.957 



^ Without IN band system. 



applying Banner reproduce the measured values from Nagano et al. very 
accurately and the partly differing deactivation constants from Bunner and 
Morozov et al. do not affect the final result much. However, this holds only for 
the comparison of the whole wavelength region between 300 and 400 nm. One 
difficulty in the measurements is the treatment of interference filters which 
have a bandwidth of about 10 nm [30]. The 10 contributions of Nagano et 
al. are given after subtracting additional contributions by smaller emissions 
within one filter region. Thus, for a direct comparison, one has to take into 
account only the 10 wavelengths reported in [30] and in this case, the calcu- 
lations with s\^Bunner differ by approximately -7%. For a detailed comparison 
of each individual band system, see Fig. 2. 

The measurements from Nagano et al. show two contributions which are con- 
siderably larger than all other data for those band systems at 315.9 and 
394.3 nm. The largest relative difference occurs at the wavelength 328.5 nm 

where the value measured by Nagano et al. is higher by 512% compared to 
the preferred calculation presented here. However, the absolute contribution 
of this band emission is only of minor importance for the entire wavelength 
region between 300 and 400 nm. 

The calculations based on the Einstein coefficients given by Gilmore et al. lead 
to typically 20% - 30% lower fluorescence yield than our preferred calculation. 
Interestingly, the band emission at 328.5 nm, which is very bright in the mea- 
surements by Nagano et al., is much lower in the calculations using Gilmore 
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Fig. 1. Fluorescence yield spectra of several calculations and measurements for 
0.85 MeV electrons as exciting particles. The bars indicate the preferred calculation 
presented in this article. All calculations are labeled with "BK_name" , where name 
indicates the authors of the input parameters. Pink stars indicate the positions of 
possible contributions of the IN system beyond the IN (0-0) transition. 
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Fig. 2. Relative comparison of 19 band systems of the calculation preferred in this 
article with measurements and further calculations. The absolute fluorescence yield 
of these contributions can be seen in Fig. 1 with the same marking. 
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et al. combined with Fons et al. data than in our modeL These uncertainties 
might be caused by possible contributions from the IN system, see Fig. 1 or 
Tab. 3. 



Hirsh et al. [33] have performed measurements for the IN (0-0) band system 
of nitrogen. They found a value for the fluorescence efficiency £391.4 nm 
0.475% which is considerably higher than the value given by Bunner, see 
Table 1. Additionally, the coUisional cross section of nitrogen with nitrogen 
and nitrogen with oxygen have been investigated [33]. The values are aisf^ = 
6.5x10^^^ m^ and (Tno = 10.9x10^^^ m^. Calculating the fluorescence yield 
for a 0.85 MeV electron with these parameters, the value at sea level in the US 
Standard Atmosphere amounts to 0.377 Pii2!225. A comparison of this number 
with the entries in Tab. 4 shows that ^1391.4 nm measured by Hirsh et al. is 
larger than our calculation by 20%, larger than measurements from Nagano 
et al. by 25%, and larger than measurements published by Bunner even by 
34%. Only measurements performed by Davidson and O'Neil result in a 17% 
higher ^1391.4 nm than Hirsh et al. 



Concluding, it can be stated that the calculations shown here provide a rea- 
sonable way of describing fluorescence emission in air while allowing for vary- 
ing atmospheric conditions. This procedure can easily be implemented into 
air shower reconstruction programs. The overall agreement in the wavelength 
region between 300 and 400 nm with some measurements is already satisfy- 
ing. However detailed, spectrally resolved comparisons reveal uncertainties in 
measurements and the understanding of the processes in air. Further inves- 
tigations are necessary, as both atmospheric transmission and EAS detection 
are wavelength dependent. Also the dependence on altitude is different for the 
emission bands. Thus, for reconstructing the fluorescence emission of an EAS 
in the atmosphere from the measured photons, all these processes must be 
understood. It must be stressed that, for the EAS experiments, the uncertain- 
ties in fluorescence yield are directly linked to the uncertainties in the primary 
energy of cosmic rays. 



It should also be mentioned that further fluorescence yield measurements, 
however without spectral resolution, can be found in literature. Kakimoto 
et al. provide a formula for calculating the fluorescence yield between 300 
and 400 nm, which gives at sea level 3.275 2^^^ [34]. This value is smaUer 
by 10.3% compared to our preferred calculation. The HiRes Collaboration 
uses a value of about 5 p^°^°"^ per charged particle in an air shower [35]. 
For these charged particles, an average energy deposit of 0.22 MeV/kg m~^ 
is assumed [14], which leads to a corresponding fluorescence yield at s.l. of 
3.811 2^^2*22^ for a 0.85 MeV electron. Assuming that the HiRes value refers to 
5 km a.s.L, one would obtain 3.6 - 3.7 ^^^^ at s.l. 
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4 Dependence on Atmospheric Conditions 



4.1 Altitude Dependence 



Firstly, the altitude dependence of the fluorescence efliciency e e^^^ in units of 
photons per MeV of deposited energy, as described in Section 3, will be shown. 
Equations (7) and (8) contain the altitude dependent variables p and T. Com- 
bined with the deactivation constants, Tab. 2, a different altitude dependence 
for the IN (391.4 nm) and 2P system (all other wavelengths) is expected, 
which is visualized in Fig. 3. Within the 2P system of nitrogen, there are no 




Fig. 3. Fluorescence efliciency profiles for different wavelengths in the US-StdA. 

differences in the altitude dependence for all emission bands belonging to the 
same excited state u' (e.g. 337.1 nm and 357.7 nm) and only small variations 
are between bands of different u' states (compare e.g. 315.9 nm from u' = 1 
and 337.1 nm from i/' = 0). With increasing altitude, the efficiency becomes 
larger due to lower rates of coUisional quenching. This increase is largest for 
the 391.4 nm band. At sea level its contribution to the total spectrum amounts 
to 8.6%, at 20 km a.s.l. it is already 10.7%, and at 30 km a.s.l. 16.8%. 

However, regarding EAS, the rate of emitted photons per meter traversed 
matter of the EAS is the observed variable which is the fluorescence yield FYx, 
Eq. (9). Via the altitude-dependent air density, Pair, the number of excitable 
nitrogen molecules and quenching partners are considered. For simplicity, FYx 
vs. altitude is shown in Fig. 4 for a 0.85 MeV electron. The most relevant 
altitude range for EAS is between ground and about 13 km a.s.l. For EAS 
with energies of about 10^^ eV, the shower reaches its maximum between 2 and 
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Fig. 4. Fluorescence yield profiles for a 0.85 MeV electron with vertical incidence in 
the US Standard atmosphere and measured Argentine atmospheres as given in [1]. 
The given yield is a sum of all emitted photons between 300 and 400 nm calculated 
as described in Sec. 3 with e^^^^^^er Table 2Morozov 

8 km a.s.l. depending on the type of the primary particle and the inchnation 
angle of the EAS. 

E.g. for the Auger experiment, the field of view of a telescope covers an altitude 
range between 0.7 km and 12.5 km above the altitude of the Pierre Auger Ob- 
servatory, 1.4 km a.s.l., at a distance of 20 km. The fluorescence yield plotted 
in Fig. 4 is a sum of all emitted photons between 300 and 400 nm calculated 
as described in Sec. 3 with ^A.Bunner ^.nd Table 2Morozov Additionally to the 
altitude dependence, also the seasonal dependence for actual atmospheres as 
obtained at the southern site of the Pierre Auger Observatory [1] can be seen 
in Fig. 4. From ground level to altitudes around 10 km, the fluorescence yield 
increases slowly. Above 10 km, the yield decreases disclosing the sensitivity to 
temperature and pressure variations. During winter I, the lower temperatures 
compared to the other atmospheric models below 9 km a.s.l. induce a higher 
fluorescence yield. Up to 17 km, the temperatures are comparatively high 
leading to a reduced fluorescence yield. During spring, summer, and autumn 
the temperatures are higher than in the US-StdA, therefore the fluorescence 
yield is decreased mostly in summer. Above 15 km a.s.l., the very low temper- 
atures during summer result in a very high emission. The differences of FYx 
for the Argentine seasons compared to the US-StdA are well below ±5%. At 
the altitude of the Auger detectors, 1.4 km a.s.l., the increase in fluorescence 
yield during winter I is negligible, however the decrease in summer amounts 
to 2.9%. At ~ 8.5 km, the differences of summer and winter I to the US-StdA 
are of the same size but with opposite signs. In winter I, FY\ is 1.5% higher 
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than in the US-StdA, and in summer 2.2% lower. More than +4% difference 
from Argentine summer to the US-StdA emerges above 16.5 km a.s.L Simi- 
lar seasonal variations in FY\ are also valid for other EAS experiments since 
similar atmospheric conditions have been found at different places [1,36]. 

The calculated altitude dependence can be compared with parameterizations 
given by authors from fluorescence emission experiments. The functional forms 
of these parametrizations are inspired by the same equations as introduced in 
Sec. 2.1 [30,34]: 



FY. 



[30] 



FY. 



[34] 



300-400 nm 



dE 
dX 
dE 
dX 



A, 



A, 



1+ pBiVT^ 1 + PB2VT 



(10) 



:iii 



While Nagano et al. [30] list A and B parameters for each of their 10 wave- 
lengths between 300 and 400 nm, Kakimoto et al. [34] just provide one set 
of parameters Ai^2 and -81,2 for the total fluorescence yield between 300 and 
400 nm. Both approaches predict similar height dependences, see Fig. 5. To 
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Fig. 5. Fluorescence yield profiles for a 0.85 IVIeV electron in the US-StdA. Com- 
parison of the altitude dependence calculated by the described method with two 
further parameterizations. 

work out the difference due to the altitude dependence, the profiles can be 
shifted so that all curves start with the same value at sea level. Then the 
parameterization by Nagano et al. agrees very well with the calculation intro- 
duced in this paper. Up to 14 km, the discrepancy is below 1%, increasing up 
to 2.9% at 20 km a.s.l. The simplified parameterization given by Kakimoto et 
al. disagrees already above 6.5 km with the calculations shown here by more 
than 1%. The difference increases up to 4% at 20 km a.s.l. 
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4-2 Humidity Dependence 



All calculations and measurements shown above are based on dry air condi- 
tions. However, in actual atmospheric conditions, there is sometimes a con- 
siderable fraction of water vapor. Thus, the effect of quenching due to water 
vapor has to be investigated. Fluorescence emission by water vapor is not 
expected. 

Our first calculations are based on Eq. (8) in which an additional term is in- 
serted to account for the collisions between nitrogen and water vapor molecules. 
The experimental determination of collisional cross sections between nitrogen 
and water vapor, which are needed in that equation, is very difficult. Two ex- 
periments have recently begun to investigate the effect of water vapor [29,37]. 

The effect of quenching due to water vapor has been studied in our preferred 
calculation for the 337.1 nm emission band. Applying the parameters from 
Tab. 2 by Morozov et al. and assuming 100% relative humidity, the emission 
at sea level is reduced by approximately 20%, at 4 km a.s.l. by roughly 5%, and 
at 8 km a.s.l. just by 0.3%. Since fluorescence telescopes typically operate only 
during "good weather" periods, this decrease in fluorescence yield should be 
considered as an upper limit. For realistic atmospheric conditions, a reduction 
of about 5 to 10% near ground and 1 - 3% at 4 km a.s.l. can be expected. 



5 Summciry and Conclusion 

EAS experiments applying the fluorescence technique measure the light emis- 
sion in air induced by charged particles, mainly electrons and positrons. The 
detected light track is converted into a longitudinal shower profile and finally 
to the total energy of the primary particle of the EAS. Therefore, the fluores- 
cence light yield has to be known precisely including spectral resolution and 
dependent of atmospheric conditions. 

The results on fluorescence yield which can be found in literature differ con- 
siderably. Most important are the fluorescence efficiency of the contributing 
band systems of nitrogen, but also the radiative life times and the collisional 
cross sections of nitrogen with nitrogen and nitrogen with oxygen have to be 
known. Up to now, a thorough understanding of the energy-dependent excita- 
tion processes of the different nitrogen states is missing. First studies can be 
found in Blanco and Arqueros [38] . 

In this article, an atmosphere-dependent description of the fluorescence light 
emission in air has been presented. The different contributions of the 2P and 
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IN band systems of nitrogen have been calculated in detail. The calculations 
are based on several parameter sets and have been compared with fluorescence 
yield measurements performed by several authors. The calculations reproduce 
some results of measurements well, while other data are off by more than 
10% regarding the total yield between 300 and 400 nm. The differences for 
individual emission bands are much larger. 

The variation of the fluorescence yield with changing atmospheric conditions 
has only been studied by a few authors. Generally, it is assumed that the main 
reduction of light emission is due to coUisional quenching. The calculation 
of altitude-dependent profiles of FYx presented here agree within 4% with 
parameterizations of measurements. 

Using the calculation and parameter set preferred in the article, a prediction 
of the influence of water vapor has been made. For realistic atmospheric con- 
ditions, an effect of about 5 to 10% near ground and less than 3% at altitudes 
around 4 km a.s.l. can be expected. Only lately experimental studies of the 
quenching rate of water vapor have been begun. 
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